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The stability of a non-ohmic/ohmic fluid interface in the presence of a constant gravitational field 
and stressed by a vertical stationary electric field with unipolar injection is studied, focusing on 
the destabilising action of the electric pressure when charge relaxation effects can be ignored. We 
use a hydraulic model, whose static equilibrium condition is written and analysed as a function of 
the ohmic fluid conductivity when subjected to a non-linear perturbation. The combined action of 
O |' the polarization and free interfacial charges on the pressure instability mechanism is also analysed. 

D . The results show some important peculiarities of the fluid interface behaviour in the presence of a 

. stationary space charge distribution generated by unipolar injection in the non-ohmic fluid. 
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I. INTRODUCTION 



On 

i 

If a stationary electric field E, parallel to a constant gravitational field g, is applied on a system composed by two 
immiscible fluids with different mass densities, the interface between them should rest, in the state of equilibrium, 
. — , , completely plane and perpendicular to the fields and eventually subjected to the destabilization when the correspond- 
ing electric field is strong enough. In general, in this paper, the term non-ohmic/ohmic, for example, refers to an 
f-H , interface where the lower fluid layer is ohmic. 

Taylor and McEwan [lj studied the static equilibrium of this system in the case of a non-conducting/conducting 
ohmic interface and determined the instability criterion in a linear theory. In this case, the stress acting on the 
interface is acting only in its normal direction: we say that the instability is due to a pressure mechanism. Melcher & 
Smith Q studied the stability of an ohmic/ohmic interface stressed by a vertical electric field in a more general linear 
theory considering all possible conductivity values of both fluids (and also other physical properties involved, such as 
^C) ' viscosity, etc.), including charge relaxation effects under these conditions. In this work, shear stresses are involved 
and hence overstability Q and surface charge convection may occur. We say that in this case the instability is due to 
ON ■ the convective mechanism. (Please note that in this case convection is due to surface charge not to volume charge, 
\ like in the electrohydrodynamic instability due to unipolar injection in an insulating liquid layer 0,13)- 

Some recent work on this problem has investigated a non-ohmic/ohmic fluid interface when the non-ohmic layer 
is subjected to unipolar injection 0, Q- These works are motivated because in certain experimental systems an 
electrode may act not only as a voltage source (a surface charge source, in the end) but also as a space charge source 
•<— 1 ' 0- Melcher & Schwartz [j| noticed that an electrode may cause, if in contact with a very low conducting fluid, 
dielectrical breakage and generate a stationary electric field with a space charge distribution in the non-ohmic part 
^ (— I of the system. This makes the coupling of the electric field with mechanical fluid equations very different from that 
occurring in the classical studies of ohmic/ohmic fluid interfaces. A clear example is the experiment by Koulova- 
Nenova, Malraison & Atten where a moderate injection in a liquid mixture of ciclohexane with TiAP salt [T(ij ] 
was produced. They observed that unipolar injection from the electrode may produce not only convection in the bulk 
of the fluid but also an interfacial instability similar to that occurring in the absence of space charge [| but with a 
peculiarity: the voltage thresholds for instability are systematically reduced by 1/3. The complete linear theory for 
a non-ohmic/ohmic interface is presented in the previous work by Vega & Perez p|, where a transition region in the 
critical behaviour of the interface has been found. This region marks the conducting-to-insulating transition in the 
behaviour of the non-ohmic/ohmic interface. In fact, the existence of this transition region implies that the dynamics 
of the same non-ohmic/ohmic interface subjected to unipolar injection may behave like in an ohmic/ohmic interface 
in which the lower layer is the most conducting, but also like in an ohmic/ohmic interface in which the lower layer is 
the least conducting (insulating behaviour). However, by definition, in the non-ohmic/ohmic interface the lower layer 
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is always the most conducting. The apparent contradiction comes from the fact that, under unipolar injection, the 
electric conduction in the non-ohmic layer may be actually more "effective" than the ohmic conduction in the lower 
layer, depending on the value of the applied electric potential. This causes the mechanics of the fluid interface to be 
much more complex when there is injection. 

The purpose of this paper is to demonstrate that this complexity appears already in the electric pressure instability 
mechanism; i.e., the static interfacial equilibrium between electric and gravitational forces. In order to make clear 
an intuitive visualization of this equilibrium, a hydraulic model is developed (figure The system described in the 
model is not that of an infinite fluid interface and thus the results are not, in general, quantitatively applicable to 
the infinite interface problem. However, as it will be shown, the results of the hydraulic model can account for the 
same transition region described above 0- Thus, the hydraulic model yields a qualitatively analogous description of 
the corresponding problem in the infinite interface. But the results are not only restricted to comparison of results 
in a previous work (7|. They also provide the following additional inputs: a) identification of the mechanism that 
causes the transition region to appear (electric pressure due to polarization charges) and precise determination of 
the transition region, and b) observation of this transition as a function of the perturbation amplitude (not only as 
a function of conductivity as previously detected Q). As we will see this implies that, once the instability begins, 
an interface with a very conducting lower layer may be stabilized in states with non-zero perturbation amplitudes. 
This behaviour differs from that observed in the purely ohmic case 0, Q , where the interfacial perturbation grows 
continually because the electrical pressure mechanism is always actively pulling up or pushing down the interface. 

The interest of this problem is due to the original behaviour of this putative fluid interface and its possible industrial 
applications. For example, the formation of stable metallic liquid points when the interface changes from conducting 
(the electric pressure has opposite sign to the applied field) to non-conducting behaviour (the electric pressure keeps 
the same sign that the applied field) in a perturbed state could be used to make ion sources. 

The possibility of producing ion sources by manipulating a fluid interface with electric fields motivated the work by 
Neron de Surgy , who extended the original work by Taylor and McEwan [| introducing a non- linear perturbation 
in a non-conducting/conducting ohmic fluid interface, but always without injection. The results proved theoretically 
that a metallic liquid can never develop stable points, independently of the geometry of the system. However, in some 
rare cases he experimentally observed stable metallic points, which Neron indicated could be due to impurities in the 
liquid. We suggest in this paper that this is related to an injection from the metallic liquid points to the air. 

In till Al we describe the hydraulic model and find its non- linear stability condition. In Mil HI we find the difference 
between the stability condition in the hydraulic model and the one found for the infinite interface in a previous 
work 0. The effect of the combined action of polarization and free surface charges will be explained in Mil CI In 
Mil Dl we introduce the reduced critical parameter Unl and the possible general behaviours of the non- linear critical 
curves are described. Finally, in MIIII and MIVI we present the results of the hydraulic model in the cases of an 
ohmic/ohmic interface and the non-ohmic/ohmic interface, respectively. Although the purely ohmic interface has 
been studied extensively 0,El^2> ^he results of MHII are interesting to make evident the new perspective gained with 
the hydraulic model. 

II. HYDRAULIC MODEL 
A. The system and the equilibrium equation 

The system (figure consists of two identical rigid cylinders with parallel axes, connected to each other through a 
thin horizontal cylindrical pipe at their base. The system is in a constant gravitational field with a constant acceleration 
g = gu z (this field acts in a direction called "vertical"; thus, its perpendicular plane defines the horizontal directions). 
The system is immersed in a fluid u, which supplies a constant pressure on another fluid (we call it fluid /), which 
is denser than fluid u. Both fluids are immiscible and incompressible, so in the equilibrium state the fluid I layer is 
below the other one. The radius r of the horizontal thin pipe is large enough to allow a negligible Poiscuillc effect for 
any typical fluid velocity (i.e., Ap/ (8filv/r 2 ) <C 1, where n is the dynamic viscosity, I is the pipe length, v the fluid 
velocity and Ap is a typical pressure difference in the vertical direction) and thus the pressure from a vertical cylinder 
is completely communicated to the other one. There is a pair of horizontal rigid electrodes in each cylinder, one at 
the top of the cylinder (at z — —L), and the other at its base (at z — d). The length of the system (L + d) is much 
lower than the horizontal dimensions of the cylinders so the boundary effects are negligible. The upper electrodes 
are connected to the same DC voltage source, that supplies a voltage V, while the lower electrodes are grounded. 
Additionally, if fluid u is non-ohmic, a space charge source at the upper electrodes injects unipolar volume charge go 
at z = —L. 
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FIG. 1: The system studied. A plane perturbation of amplitude rj is introduced in such a way that the interface level is raised 
in the cylinder centered in x = and lowered in the other cylinder (which is centered at x — A/2). 



We write the Navier-Stokes equation: 

= -Vj> + /iV 2 v + V • T e + pg (1) 

where p is the mass density, v the fluid velocity, p the total pressure and T e denotes the electric stress tensor, whose 
elements are 7% = eEiEj — ^SijeE 2 (s is the dielectric constant of the fluid, subscripts i and j indicate the components 
in cartesian coordinates (i,j=x, y, z) of the stress tensor and the electric field of modulus E, and 5^ are the elements 
of the identity tensor) . 

We define the jump of a magnitude A as the difference between its values just below and just above the interface 
and denote it as < A >= Ai(F(r)) — A u (F(r)), where F(r) is the interface position. The normal stress balance 
condition in the interface is written: 



n • (T v + T e ) n - (p) = 7 V • n (2) 

where n is the normal direction to the interface, 7 is the coefficient of surface tension and T v is the viscous stress tensor 
for incompressible fluids: = \x(dvijdj + dvj/di) (again, i 1 j=x 1 y, z). From now on, we will ignore coordinate y, 
due to the system symmetry. 

If the system is in equilibrium, the fluid interface is horizontal and at the same level in both cylinders (figure^. 
In this situation all stresses at the interface are in the vertical direction and we have v = and V • n = V • u 2 = 0. 
Then the normal stress balance reads: 



(P) = (^ 2 ) (3) 

We introduce now a plane perturbation and the interface level is raised a height rj in one cylinder (the one at 
x = 0) and lowered the same height rj in the other (Fig. 0. The perturbation is kept by some pressure source until 
the electric field becomes stationary. In this point, we still have v = and V ■ n = V ■ u 2 =0. Now the pressure 
source stops and the system is left under the action of gravitational and electric pressures. It is implicit in the way 
in which the perturbation is introduced that at this point charge relaxation may be ignored and that condition @ 
is still fulfilled. We express now the pressure as a function of the scalar field II, which is defined by the relation: 



II = p — pg ■ r and that is called " modified pressure" [13j . The equation in an equilibrium state may be written 
then: 



vn = V • T e (4) 

which reflects the fact that the total pressure that II the surface a body immersed in the fluid feels is p modified by 
the gravitational force 0] . When the gradient of this modified pressure II is in balance with electric stresses in the 
volume of the fluid, like in flf. there is no net volume force in our system. Respect to the net pressure jump at the 
interface, it can be rewritten as a function of the modified pressure IT: 



(pL-(p)OT=(n)_, {p) + -(p)g(- V ) = (U) + (5) 

where subscript "— " stands for the value of the magnitude at the interface in the cylinder at x = A/2 (downward 
perturbation) and subscript "+" stands for the value of the magnitude at the interface in x = (upward perturbation). 
If we use these expressions into Q and we take the difference between the pressure jumps in both cylinders we obtain: 



II II - • {If^j -(fE 2 ) 

When the modified pressure jump at the interface in both cylinders is the same, the surface forces at the interface 
in both cylinders are equilibrated. Thus, our stability condition is: 



(n)_-(n} + = o (7) 

This is the only mechanical equation we are actually using in this work, from now on. When the electric pressure 
term overcomes the gravitational term then (II) _ — (II) , > and the perturbation tends to amplify. If (II) _ — (ft) + < 
the perturbation tends to damp. The stability condition leads Q) to: 



(l £E2 ) -(^ E2 ) + = 2 (P)9V (8) 

In the absence of electric forces, the stability condition JSJ gives the solution 77 = (i.e., if no electric field is applied, 
logically the only possible equilibrium state is the same interface position in both cylinders). 

It can be demonstrated that for infinitesimal 77, the hydraulic model stability condition JSJ reproduces the linear 
instability criterion for an infinite plane fluid interface under a vertical electric field in the limit of long wavelength if 
the horizontal variation of the pressure is zero, as we will see in the following section. 



B. The hydraulic model and the problem of the linear stability of an infinite plane fluid interface 

In the analogous and more studied problem of an infinite plane fluid interface stressed by a vertical stationary electric 
field, several mechanisms may simultaneously induce an electrohydrodynamic instability, unlike in the hydraulic model 
here developed, where the electric pressure is the only destabilizing mechanism. In order to make a comparison between 
the results that will be presented in this work and the results in the problem of the infinite interface, two questions 
could be posed: 1) are there situations in which, like in the hydraulic model, the electric pressure is the only active 
instability mechanism in an infinite plane interface?; and 2) if so, can the hydrostatic model account for its linear 
instability threshold values? 

Concerning the first question 7] demonstrated that in the infinite plane interface the electric pressure is the only 
destabilizing mechanism involved when an initial perturbation with infinite wavelength (A = 00) was produced. And 
this infinite wavelength instability occurs if capillary forces are strong enough 1] and additionally, in the specific case 
of a non-ohmic/ohmic interface, if the non-ohmic fluid has a very high ionic mobility Q. Assuming that we deal 
with an infinite plane fluid interface whose properties fulfill these conditions (i.e., that the long wavelength instability 
occurs) and concerning question 2), the answer is yes, although not always. As we will see, this is because in the 
hydrostatic model the variation of the pressure in the horizontal is not taken into account. This can be shown if 



the linear perturbation of Navier-Stokes equation term for horizontal components x, y are considered. Given the 
symmetry of the system, it is enough to analyse the x component: 



d5p t 35(f) ( 2 ® 



where q is the free space charge in the equilibrium state, and 5<p and Sp are the electric potential and the total 
pressure linear perturbations in a very slightly deformed interface, while Sv x is the x component of the velocity linear 
perturbation. 

Integrating this equation in x and taking into account that Sv x — in the limit of long wavelength .7] , when the 
jump at the interface is taken the following relation is obtained: 

(SU) = - (qSct>) (10) 
Consequently, the total jump of the modified pressure at the interface is now: 

(II + SU) = (^£ 2 ) - (P) 9*. - (qS</>) (11) 

where z s is the interface level. Thus 

The stability condition to be used for an infinite interface is the following: 



(n + (5n)_ - (n + (5n) + = o (12) 

In the case of an ohmic/ohmic interface the additional term is always zero because no volume charges are present 
(q — 0), which means that the hydraulic model yields the exact linear criterion for an infinite plane interface in the 
limit of long wavelength. However, in the non-ohmic/ohmic interface this term is not zero except in the cases of a 
perfect conducting ohmic fluid, ai — oo (being 07 its electric conductivity), and a non-conducting ohmic fluid, ai = 0, 
where the interface is an equipotential, and thus 8(f) — in the interface. It is convenient to comment at this point 
that (|12fl is equivalent to the linear stability condition for the infinite interface found in an author's previous work 
0, as we will numerically check out in section 4. The advantage of starting out from the hydraulic model is evident 
if we notice that the deduction of the stability condition has become now much simpler ■ 

As we see, then, the additional term (SU) is in general needed to obtain the exact criterion for an infinite and 
initially plane non-ohmic/ohmic interface. Nevertheless, the simpler stability condition J7J obtained for the system of 
two cylinders not only describes essentially, in a non-ohmic/ohmic interface under unipolar injection, the instability 
regions as a function of the ohmic conductivity but also the linear criterion threshold value, as we will see in £IIVI 

We recall that the comparison to the hydraulic model is restricted to an infinite wavelength perturbation in the 
infinite interface. For shorter characteristic wavelengths capillary forces are involved but the present analysis is useful 
because it provides an intuitive description on the mechanical process that occurs at the interface due to the action of 
electric pressure against gravitation. And this action is present as the main destabilizing mechanism in any problem 
of a fluid interface stressed by a vertical electric field. 



C. Dimensionless magnitudes. An intuitive framework 



We take d, (p) gd, and \J (p) gd 3 /e u respectively as reference units for distance, pressure, and electric potential, 
being e„ the permittivity of fluid u. From now on we will only use non-dimensional magnitudes, and we will denote 
them with the same symbols that we used for the dimensional ones, except for the perturbation amplitude, that we 
will call £ = -q/d. 

The static equilibrium of the interface in the hydraulic model is given by the opposition of a gravitational term and 
an electric term. In a perturbed state the gravitational term always acts towards the part of the system with a lower 
thickness of the heavier fluid layer. On the contrary, the electric pressure may act towards any of the two cylinders, 
depending on both the magnitude and sign of the electric pressure jump in each cylinder. Thus, it is convenient to 
write the dimensionless electric pressure jump in the interface using a parameter that allows to easily determine case 
by case the electric pressure sign: 

(¥ E2 ) = h El{Zs)2 - \ E ^ z °? = \ E &°¥ & - s2 ) ( 13 ) 
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FIG. 2: A system with E < Eo when the pressure instability mechanism is possible (E < y/el). Long white arrows indicate 
the action of electrical pressure jump. In (6), when a perturbation is introduced, we see a gravitational force appear. This 
gravitational pressure flow, directed towards the cylinder with lower interface height in all cases, is 2 (p) grj (or 2£ in reduced 
magnitudes). In this case, the net electric pressure flow is directed towards the left cylinder (three arrows against one). 



where we use subscripts u and I to denote the magnitudes in fluids u and I respectively, z s is the interface position 
and £ = E u (z s ) / Ei(z s ) > 0, that we call the "apparent conductivity" of the interface. It reflects which of the two 
fluids is more conducting in the interface: £ > 1 if E u (z s ) > Ei(z s ) and we say that the interface is conducting and 
conversely, not conducting if E < 1. The expression is valid for all fluids independently of their regime of electric 
conduction. 

The total surface charge Qt an d the free surface charge Q at the interface are, respectively: 



Qt=e (l-H)E l (z a ) Q=(e l -E)E l (z s ) (14) 
where eo is the reduced vacuum permittivity. The stability condition © in reduced magnitudes gets: 

(15) 

There are two special cases for which an eventual instability due to electric pressure is not possible, as the left hand 
term in i|15fl is null. In effect, in the case £ = Je[ it is evident, from l|13fl . that each one of the electrical pressure 
terms in the left hand of Ijl5(l is zero so there is no non-trivial solution (£ ^ 0) to this equation. And in case the 
equality (e_E 2 )_ = (e_E 2 ) + is fulfilled, the left term of the stability condition l|15fl is again null and a non-trivial 
solution does not exist. 

The electric pressure jump, or equivalently (|13fl . Ei(z s ), decreases with the lower layer thickness if the apparent 
conductivity is high enough; i.e., the lower layer is "less conducting". This occurs if E < Eo, where the value of Eo 
depends on the electrical regime of conduction of the fluids. Figure[2Ia) represents the initial unperturbed state when 
the electrical pressure jump is positive (E < ^/ei) and decreases with the lower layer thickness (E < Eo). 

If a plane perturbation is introduced, the magnitude E gets in general a value E_ in the cylinder with a minimum 
interface elevation and another value £+ in the cylinder with maximum elevation. But in order to simplify, let us 
restrict to an infinitesimal perturbation, so the new values still fulfill E± < So 5 Once the perturbation is introduced, 
a gravitational pressure acts against it by communicating an upwards pressure to the zone with minimum interface 
elevation (figure^). Now the electric pressure jump (|Ll|l is higher in the zone with a lower interface height as the 
thickness of the lower fluid layer has decreased. Conversely, in the zone of maximum elevation the electric pressure 
jump decreases. Thus, a net electric pressure flow appears towards the left cylinder in figure b). If the difference 



a) 



b) 



d+r) 



f t 


i g 


t f 






i 


1 


+ + + + + + + + + + 

t 


jg t g 


+ + + + + + + + + + 

Itt 







I 



r 



d-n 



FIG. 3: A system with E < Eo when the pressure instability mechanism is not possible (E > y/ei). We see now the net electric 
pressure flow is directed towards the right cylinder (lower interface height), in the same restoring direction that the restoring 
gravitational term: the electric pressure is stabilising. 

between the electrical pressure jump in both cylinders is high enough, the restoring action of the gravitational pressure 
will be counterbalanced. This is possible if the applied potential is higher than a critical value V c , given by the stability 
condition Q. 

On the contrary (figure UJj), if the electrical pressure jump magnitude still decreases with the lower layer thickness 
but becomes negative {^/£l < E < Eo), the perturbation is damped (figure The same analysis may be carried 
out when Ei(z s ) increases with the layer thickness of fluid I (E > Eo), so finally two linearly stable regions are found: 
y/ej < E < Eo and Eo < E < ■Je[. Thus, the apparition of linearly stable bands is related to the change in the 
tendency of the electric pressure with the lower layer thickness, which gives the limit Eo, and to the change of electric 
pressure jump, which gives the limit £i. 

We see that for a given behaviour of E[(z s ), the stabilization is due to a change in the sign of the electric pressure 
jump. In order to find out what may cause this change of sign, let us write the electric pressure jump as a function of 
polarization (Q p = Q t — Q) and free interfacial charges. The contribution of the free surface charges to the electrical 
pressure jump has the same sign of the total interfacial charge term: 



E 

—Qt + Q 



(16) 



while the contribution of the polarization charges has opposite sign to the term of total charge: 



\Ei{z a ) 



l + — )Qt~Q P 



(17) 



Thus, we see that if no polarization charges are present initially, the electric pressure jump takes the sign of the total 
interfacial charge Qt- If we now change e; at constant E in such a way that polarization charges have the same sign 
that Qt and get high enough, they can change the sign of the electric pressure jump. Thus, they play an essential 
role in the stabilization of the interface. 



D. On the critical curves in the hydraulic model 



It is also convenient to define the dimensionless parameter Unl — £uV r c 2 /((p) gd 3 ) where V c is the applied potential 
for which the stability condition Q is fulfilled (i.e., V > V c yields (IT) _ — (11), > and the perturbation can be 



sustained). The parameter Unl represents the critical electric pressure e u V 2 /d 2 , reduced with the gravitational 

1/2 

pressure (p) gd, while its square root U^ L represents the reduced critical applied potential. Un l is in general a 
function of the perturbation amplitude £. 

Let us start with an unperturbed state (£ = 0, V 2 < Unl{£, =0)), where an arbitrary stationary perturbation with 
amplitude £o is introduced. Then a valid solution Unl((,o) from (|15|l should fulfill two conditions to make possible 
the perturbation be sustained: 

*; u NL (Zo) > o. 

ii)V 2 >U NL (Z ). 

The first one comes from the definition of Unl because [//vx(£o) < should correspond to an imaginary critical 
potential. The second one is also necessary because V 2 < f/jvx(£o) indicates that the perturbation is decreased to a 
lower value. 

Typical critical curves Unl(0 are represented in figures a) andJ^J In these curves stable and unstable regions 
are delimited by the function Unl(£,) and the behaviour of Unl(£,) provides information about the evolution of the 
perturbation once it is introduced. An increasing Unl{C) in £o {dU(£ t o)/d£ l > 0) means that the perturbation will not 

tend to increase for an applied potential V = ^vl(£o), because for £ > £o> V < ?^vx(0 (the interface cannot overcome 
the restoring action of the gravitational pressure). And viceversa, a decreasing Unl in £o means that the perturbation 

will tend to increase if V = ^vx(£o)- The third case occurs when d[/(£o)/<9£ = 0. In this case, for V = E/jvl(£o) the 
perturbation will not increase if d 2 UNL{£o)/d£ l 2 > but it will tend to increase if d 2 V 'n l(£,o) / 'd^ 2 < 0. In definitive, 
if dU NL {S,o)/dS, > (or if <9[/(£ )/<9£ = 0, d 2 U N x(£o) / 'd? > 0) it may be said that "stabilization" of the perturbation 

1/2 

occurs at V — U^^q) and the point £ = £o may be considered as a new point of stable interface position, that is 
different to the trivial solution £ = of the initial equilibrium state. We call these points "perturbed stable states" 
and as we will see they are only present in the non-ohmic/ohmic interface. 



III. OHMIC/OHMIC INTERFACE 

In the case of two ohmic fluids with conductivities a u , er; , we take a u as unit for electric conductivity. Taking into 
account that in the stationary state V • j = 0, then j u = ji = j, E u = aiEi and the expressions for the dimensionless 
stationary electric field and E in the unperturbed state are the following: 

Ei = Vi E u = V u /L = aiVi/L E = 07 (18) 

being Vi and V u the electric potential drop through the lower and upper fluid layers inside the cylinders. As the 
perturbation is plane, the fields in the perturbed state are: 

E l± = E u± = (19) 

where the upper signs stand for the magnitudes evaluated at x — and the lower signs stand for the magnitudes at 
3; = A/2. 



Using V = V u + + Vi + = V u - + Vi-, we obtain 



^E") = 1 -(e l -a 2 )-^-, (20) 



Vl± = ML?\na±t) v (21) 

In an ohmic/ohmic interface E is a constant (|18(l and therefore we can make the interface at constant E coincide 
with the real interface (constant conductivities of the fluids). It is also easy to find that E = 1: E[(z s ) always 
decreases with the I layer thickness for E < 1 and always increases for E > 1 Ijl9|l . Thus, given the conductivities and 
permittivities of the fluids, the interface should be stable respect to the pressure instability mechanism for any value 
of the electric field, if ^fel <E<lorl<E< ^/ej, following the analysis in £111 CI This can be confirmed by finding 
the function Unl- 

After some short calculations and using the stability condition in reduced magnitudes l|15|l we get Unl as a function 
of the perturbation amplitude and the other parameters of the system: 



U N L((,,cri,£i,L) = f((Ti,ei,L)JC(^,a-i,L) 



(22) 
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FIG. 4: (a) The parameter Unl (if positive) is always decreasing with the perturbation amplitude; a; = 100, £; = 4, L = 1. 
(6) The linear criterion (Unl(^ = 0)) plotted against <t; presents a stable region (the negative values of Unl), as a consequence 
of the action of polarization charges; e; = 4, L = 1. 

where /(07, £7, L) and /C(£, 07, L) (which gives all the dependence on the perturbation amplitude £) are the following 
functions: 



f(<ri,£i,L) 



{ei - af)(l - a t )(l + Lat) 



1C& a u L)= ((1 + La,) 2 - (1 - aj) 2 ^) 2 (23) 

It is to be noticed that for S = ct, = 1 and S = o\ = the parameter C/jvl takes an infinite value for any £ (i.e., 
the instability is not possible) , which is in accordance with the analysis in £111 CI (0/ = ^fei is equivalent to S = -^/eT 
and 07 = 1 gives (eE 2 ) _ = (eE 2 ) + ). 

In an ohmic/ohmic interface, condition i/ in £111 Dl is fulfilled if and only if f(ai,ei,L) > 0. This is so because the 
function /C(£, 07, L) is always positive And the sign of /(cr;, £/, L) is positive l|2*3l unless the terms (1 — af/si) and 
(1 — 0;) have different signs. This occurs when ^fel < 07 < 1 or 1 < 07 < i.e., no perturbation can be sustained if 
the electrical pressure jump and the total surface charge Q t at the interface have opposite signs (see ifUjl. (|TB|) . J2DD- 
Thus, the result advanced by the analysis in £111 CI is confirmed by the explicit calculation of Unl- Analysing further 
the conditions for stabilization we notice that they are fulfilled when the interfacial free charge Q = (e/ — T,)Ei(z s ) 
takes opposite sign to the total interfacial charge Q t , and consequently the polarization charge Q p = Qt — Q takes 
the sign of Q t , which agrees again with the analysis performed in i jll CI Although the problem of the stability of an 
ohmic/ohmic interface has been extensively studied, this stabilizing effect due to polarization charges has not been 
formerly detected. The values of Unl are plotted in figure 6) vs - the reduced conductivity 07. The stable region 
corresponds to the interval where Unl takes negative values. 

Besides, we saw in £IIIDI too that for V 2 = Unl(^o) the perturbation amplitude tends to grow up to £1 > £0 from 
its initial value £0 if V 2 > Unl(Q inside the interval (£i,£o), which is always the case in an ohmic/ohmic interface if 
conditions i) and ii) are fulfilled. In effect, let us study the derivative dUNL/d^: 



dU NL TS dK{£,ai,L) 

on 



9£ 

= -4£( ffi - l) 2 ((1 + La t ) 2 - (0/ - l) 2 e 



(24) 



As we see, if f(ai,ei,L) > 0, 9C/jvl/9£ is always negative (figure^), provided that £ < L (a perturbation with 
£ > L forces the interface to touch the upper electrode, case that we do not analyse). The result is independent of 
the initial of the values of £0 and £1. This means that in an ohmic/ohmic fluid interface the pressure mechanism is 
always self-fed as it becomes increasingly stronger: once the instability is set on the perturbation amplitude grows up 
to its maximum value. This is a well known characteristic of EHD instabilities in plane interfaces in ohmic systems 

El 

For 0/ — > 00, it is easy to find from 11' 21) and l|23|l that Unl{£, =0) = L 3 , wich agrees with the result by And 
finally, the author has checked that for the linear instability at finite conductivities the critical values provided by 



the hydraulic model (Unl{£ — 0)) coincide exactly with those of the complete linear theory for an infinite plane 
interface Q in the case of an infinite wavelength instability with negligible charge relaxation effects, confirming the 
demonstration in 311 Bl 

We see then that although the instability in an infinite ohmic/ohmic interface has been extensively studied, the 
hydraulic model is able to reproduce some former basic results 0, UJ$ and also to find a new feature: the stabilizing 
effect of the polarization interfacial charges when combined with the action of free interfacial charges. But the most 
relevant new results of the hydraulic model are found in the non-ohmic/ohmic interface under unipolar injection, in 
the next section. 

IV. NON-OHMIC/OHMIC INTERFACE 
A. Equations 

Let us suppose now that the fluid " u" is in non-ohmic regime of electric conduction and the fluid " 1" is in ohmic 
regime. If there is a unipolar space charge source in the upper electrode, a unipolar space charge distribution is 
induced in the non-ohmic fluid, in which the conduction (in dimensional magnitudes) is expressed by j„ = q u (K u ~E u + 
v) — D\Jq u ] where Kyis the ion mobility of non-ohmic fluid and D u its diffusion coefficient. The diffusion term may 
usually be neglected Q, so in a state at rest (v = 0) we have in our system that j u = K u q u E u in the non-ohmic 
fluid. We take now K U ^J e u (p) g/d as unit for electric conductivity, being K u the ion mobility in fluid u. In reduced 
magnitudes we have the following electric equations in stationary regime (for which V • j = 0): 

ju = q u E u = j = q u (25) 

3 i = aiEi =3 ^ = (26) 
dz 

with the boundary conditions in the electrodes: 

M-L) = l q u (-L) = C <M1) = (27) 

and in the interface: 

0„(o) = ^(o) i«(o)=ii(o)=i (28) 

where the condition q u (—L) = C denotes the fact that there is a space charge source at z = —L. The parameter C 
represents the reduced space charge that the upper electrode injects on and is usually called "level of injection". The 
non-ohmic conduction due to unipolar injection and the correct boundary condition for a unipolar injection source 
q u (—L) have been studied in detail by Atten 0,0, and references therein for more details on this issue]. 
Then the solution of the stationary electric field is: 



E u (z) = ^2j(z + b) Ei{z) = l- E = J—a t (29) 

oi 

with b = 2^2- + L. 

We see that now E is not a constant but a function of the electric field (through the current density) and it is 
always possible for any real interface (i.e., any given value of <ji) in the initially unperturbed state to find ranges of 
the electric field for which its corresponding E lies out of the stable intervals: there exist real Unl((, = 0) which 
are solution of the stability condition J7J for all values of ui . This means that the linearly stable region around the 
intervals yje[ < E < Eo or Eq < E < ^fel only exist as a consequence of dimensionalization and therefore they should 
disappear for a real interface (constant ai ) . In any case we will see in fllVDI that now the action of polarization 
charges affects the pressure instability mechanism in other ways. 




If the system is under strong injection conditions (i.e., C — > oo), the expressions for the electric field and £ in the 
unperturbed interface are the following: 



E t (0) = Vi E u (0) = x/2jl = y/2aiViL 



2mL 
Vi 



(30) 



If a plane perturbation of amplitude £ is introduced, the equation for the electric pressure jump in the cylinder at 
the interface is: 



1 9 



V 2 
v l± 



(i±0 2 



2Vi±ai 



Lrj. 



(31) 



The electric pressure jump can be expressed as a function of the applied potential V taking into account: j u + = ji+, 
j u - = ji- and V = Vu— + Vi- — V u + + Vi+ for the electric potential. Then, in the perturbed interface we get: 



(V - V l± ) 2 = p±v h 



(32) 



where 0± = a x {L T 07(1 ± 0- 

From l{3"2l , we get the solution of Vi± : 



(33) 



The other root of (|32[1 . the one corresponding to the sign + before the square root, is not possible as it gives 
solution of Vi± such that Vi± > V. In effect, as y/ j3 2 + A(3V > (3: 



V, 



l± 



2V + P± + yJfc + W±V 2V + (3 : 



> 



> V 



(34) 



The value of So is not constant in the non-ohmic/ohmic interface (from l|31|) and l|33[) ) and in general Eo ^ 1, unlike 
in the ohmic/ohmic case. We introduce the expressions for the electric pressure jumps in the stability condition in 
reduced magnitudes l|15|) and then we get: 
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V?- 



(1-0 2 



(i-0 
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(1+0 5 



(1 + 



4£ 



(35) 



As Vi+ and Vi- are, from (|33|l . determined by V , when the equality is held in i|35|) we get the condition of minimum 
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applied potential UJ L for which the instability mechanism is possible. Unlike in the ohmic/ohmic case, now (in 
general) there is no analytical expression for Unl- 



B. Perfect conductor limit 

Let us suppose that the ohmic fluid is a perfect electric conductor (07 — > 00). In this limit V; — > and then: 



2-^7^2^ « ^± = 7^ + ^-^! + — J (36) 

where a± = (8/9) (L T 3 /(l ± 0- 

We analyse the limit of <j{V\ when u\ — ► 00. Developing the square root in power series of I/07 << 1, we get from 
and 



(37) 



being 



m - (L + 2 -(L-0 2 ~i ( e } ( } 

The function gives the variation of the non-linear criterion with the perturbation amplitude £. 

Let us compare now with the case without charge injection (C = 0). In the absence of injection the solution of the 
electric field at the interface tends to E u ± — > V u ±/ (L =p £). Operating in a similar way to the strong injection case we 
have: 



Unl = H€) (39) 

that is exactly the same to the case of strong injection (|37|l except for the factor 4/9 that now does not appear. 

The dependence of Unl (and U^f 2 ) with £ is, from (|37l 139(1 . the same for C — oo and C = and is given by the 
function (or T 1/2 {£) for U]H). We can see this behaviour of f7jvx(£) in fig. EJa): for ai — oo (and hence, also for 
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(j / =0), UJ L is always a decreasing function of £. Then, once the instability starts it tends abruptly to states with a 
minimum Unl, which are the ones having a maximum value of the perturbation amplitude. An analogous behaviour 
has been experimentally observed in a conducting liquid with Q and without injection [Hlllj|. who observed that the 
instability develops violently towards the upper electrode, producing an electric breakage. 

In the limit of zero perturbation amplitude the linear criterion for the instability in a perfect conducting fluid 
interface can be reobtained. In effect, if £ — > we obtain that — > L 3 , and then U^f 2 = L 3 / 2 for the case without 

injection, which agrees with previous works and U]J 2 = (2/3)L 3 / 2 for the case with infinite injection, which 

agrees again with the result in former works 0, ■ Note also that the case C — yields the same criterion that in 
the ohmic/ohmic case for oi — > oo, in flllll 



C. Non-conductor limit 



Analogously, the non-conductor limit can be taken. Developing l|35|) in power series of ai <^ 1, and in the limit of 
o i — ► we get: 



"»* - [7.) + = S (1 " S ' = ( " 0) 

where now before the function of £ appears a factor l/ei instead of the 4/9 in the perfect conductor. The value u; = 
has not real physical meaning but it is interesting the study of this limit as a reference for very low conductivities. In 
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definitive, the behaviour of UpJ L (t;) for ct; 1 is similar to that in the limit ai — > oo, represented in figure|3Ja)- This 
can be seen analytically in the function 7i(£), which has the same behaviour that the function Unl{£) is always 

decreasing and is zero for the maximum perturbation amplitude. A similar behaviour is detected in experiments in 
very low conducting liquids under unipolar injection: the rose-window instability has a high deformation amplitude 
(of the order of the liquid layer thickness) near the instability threshold The peculiarity respect to the high 

conducting case is that now the electric pressure is directed downwards. 

The zero perturbation amplitude limit, f — > 0, yields H(£) — 1 and U^ L = 1/y/ei, which agrees with the result in 
the linear theory for the infinite interface pfl. 



D. Non-linear transition from conducting regime to low conducting regime 

In aTVBl and CVCl we have demonstrated that in the limits of a non-conductor and perfect conductor ohmic layer 
Unl(£,) is minimum for the maximum perturbation amplitude. This means that the instability, once is set on, evolves 
up to the maximum perturbation amplitude. The difference between both limits comes from the fact that while in 
the perfect conductor the instability is due to an upward electric pressure, in the non-conductor limit the instability 
is driven by a downward electric pressure. It seems reasonable to think that between the perfect conductor and the 
non-conductor behaviour there should be intermediate behaviours, i.e., perturbations that do not grow up to the 
maximum value. In effect, we saw in £HV AI that the real non-ohmic/ohmic interface is always initially unstable: it 



is possible to find a finite value of the applied electric field that is able to sustain any finite perturbation, which is 
a consequence of the non-constancy of S l|30(l in the non-ohmic/ohmic case. But this does not prevent, after the 
instability is set on, the interface from entering a stabilizing behaviour. There would be two possibilities: one is that 
the electric pressure jump changes sign, and the other is that this pressure jump changes its behaviour with the lower 
layer thickness, while the interface is evolving. Either of the two possibilities could make the interface enter the stable 
intervals £o > £ > ^fe[ and */e7 > £ > £o. We recall this is possible in the non-ohmic/ohmic interface only because 
the apparent conductivity (|30|l is not constant, and thus the interface is taking new £ values as the perturbation 
evolves. Thus, although linear stable states do not exist now for any value of 07, it should be possible to think in these 
cases of an instability that evolves without reaching the maximum perturbation amplitude: we say that perturbed 
states are expected to appear. 

We saw in i|II Dl that the perturbed stable states occur in an interval of £ when dUMh/d^ > in it. And, in 
effect, critical points of this type arc calculated for the first time in this work, in figure \^b), where U^ L is plotted 
for a non-dimensional conductivity 07 = 10 4 . This reduced conductivity corresponds in an air/liquid interface to 
a physical conductivity of the order of 10~ 2 fi m when the dimension of the system is of about 1 cm. We see 
in figure b) that the function Unl{0 (or, equivalently Unl^) 1 ^ 2 ) obtained from has an increasing behaviour 
near the maximum amplitude. This means that a fluid much more conducting than water, for example, shows such 
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stabilization, which could be possible even for an initial value £0 = and V = U^ L (0) as once the interface reaches 

the increasing E/jy£(£) region it begins to be decelerated as soon as it takes a value £1 for which V < U]j 2 {^i). 

We saw in the linear theory for an infinite plane interface that there may exist multiple instability threshold values 
for a given conductivity [7J. This also occurs in the hydraulic model because if we further decrease the conductivity, 
three solutions at £ = begin to appear (figureJSJs). Two of them enclose an unstable region which should correspond 
with the high conducting regime-like instability (upward electric pressure jump) while the third one corresponds to 
the low conducting mechanism (downward electric pressure jump). The two solutions enclosing an unstable region 
get nearer as the conductivity is decreased (figure |5p): the conducting mechanism tends to disappear. Finally, for 
lower conductivities, a unique solution at £ = corresponds to the non-conductor limit of 1)41)0. The value of £ 13(J|) 
can be used to determine the corresponding pressure mechanism. 

Figure represents U]J 2 {^— 0) as a function of the conductivity 07. The curves from the stability condition for the 
hydraulic model (7J show the same qualitative behaviour that the linear critical values from the stability condition for 
the infinite interface 1)12(1. In addition, the linear critical values from l|12|) coincide (figureEJ) with those for an infinite 
interface in the dimensional representation (scaled to our dimensionless magnitudes) of a previous work [JJ, which 
confirms that the stability condition l|12|l applies for an infinite interface. The quantitative similarity of the results 
from (J7J and 112|) in ngureHHsuggests also that additional term in the stability condition for the infinite interface is less 
important. Thus, the limits of the stable bands detected in the dimensionless representation of the author's previous 
work must be close the limit values for the hydraulic model: y/ej and £0 Thus, the stable bands are related to 
the change of electric pressure jump sign and the change of tendency of the electric pressure jump with layers relative 
thickness, instead of being related to the change of sign in the interfacial charge, as suggested previously Q (notice, 
too, that the change of sign of interfacial charge does not provide two limits but just one value). 

V. CONCLUSIONS 

We study a hydraulic model, for which we find a non-linear stability condition, with a twofold objective: 
First, to study specifically the unstabilizing mechanism of electric pressure jump against the gravitational force, 
we use a dimensionless representation in which we introduce the fundamental parameter "apparent conductivity" 
(£). This allows a systematic study of the different types of electric pressure unstabilising mechanisms. The simple 
and intuitive hydraulic model and the use of the apparent conductivity, allow us to find two stability bands and 
their origins. This is in the action of polarization charges, that may inhibit the instability pressure mechanism in 
an interface subjected to a vertical field when combined to free interfacial charges. In our formulation, the stability 
bands are: £0 > £ > \fz~i an d ^fe[ > £ > £0, where Ei and £0 are the reduced dielectric constant and the value of the 
apparent conductivity for which the electric pressure passes from decreasing to increasing with the reduced lower layer 
thickness. As a consequence of this, the pressure instability mechanism is not possible for any ohmic/ohmic interface 
with reduced conductivity 07 in the intervals 1 > 07 > */e7 and ^fel > o\ > 1; i.e., the total surface charge and the 
electric pressure jump at the interface have opposite signs. In the case of a non-ohmic/ohmic interface, the pressure 
instability mechanism is always possible, for any value of the ohmic conductivity, when we pass to the representation 
in 07. 

In the hydraulic model there is also a transition region 7] in the pressure instability mechanism in a non- 
ohmic/ohmic interface for a linear perturbation, but also for a non-linear perturbation. The behaviour of the parameter 
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FIG. 5: Uf/l as a function of the non-linear perturbation amplitude £ in a non-ohmic/ohmic interface under strong injection 
(dimensionless L — 1, si — 4 in all figures), (a) For a non-ohmic/perfect conducting interface. (6) For a high non-dimensional 
conductivity value (ai = 10 4 ). (c) For a lower non-dimensional conductivity (cr; = 20) three solutions in £ = appear, (rf) The 
two first solutions in £ = tend to disappear if the conductivity is still lowered (07 = 7). 



Unl (or Uj L ) as a function of the perturbation amplitude reveals the existence, for intermediate conductivities, of 
perturbed stable states (which are different to the trivial solution of the initial equilibrium state £ = 0). These new 
stable states are only found in the non-ohmic/ohmic interface , su ggesting the possibility of interfacial dynamics very 
different to those detected in systems without injection [l|,@,[l2|. This could be related to the observations by 
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of stabilised metallic points. In the limits of very low or high conductivity the behaviour of the function U^ L {^) is 
as expected (analogous to the observed in the infinite interface): the instability evolves up to the maximum value of 
perturbation amplitude. 

Although a detailed study of the interfacial dynamics is needed in order to do determine precisely the situations in 
which the stabilisation in perturbed states from any initial state is possible, it seems evident that the introduction of 
the injection enriches the behaviour of the pressure equilibria in a fluid interface, and also opens a way to stabilization 
and control of the interface deformation of high conducting fluid interfaces by applying stationary electric fields. 

The second main objective of this work is to approach, in a mathematically simple way, the stability condition found 
in a previous work for an infinite non-ohmic/ohmic interface under unipolar injection, in the limit of long wavelength 
0- In this issue, an equivalent stability condition to the referred one is found in a very simple way, in HII Bl The 
difference between the linear stability conditions in the hydraulic model Q and in the infinite interface l|12|) was 
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FIG. 6: Reduced critical potential Unl vs. reduced conductivity oi for an infinitesimal perturbation £ — > with L = 1. 
Continuous curve stands for the hydraulic model stability condition Q while discontinuous curve stands for both the stability 
condition in a previous work and 11211 . (a) Reduced dielectric constant ei —4. (b) Reduced dielectric constant e; = 40. The 
effect of increasing the dielectric constant (i.e., increasing the polarization interfacial charges in the high conducting region) in 
the non-ohmic/ohmic interface is to move the conducting critical curve (the one at the right) to higher values of <t;. 

also found. This allows us to determine if there are cases in which these conditions coincide. We found that this 
coincidence occurs in the case of an ohmic/ohmic interface and also in the non-ohmic/ohmic interface in the limits 
of perfect conducting and very low conducting ohmic fluid. We also found that for intermediate conductivities, the 
difference between linear critical values given by the two stability conditions ((71 112fl is not important (figure EJ. This 
allows us to assert that the limits of the stability bands found are actually related to the values y/el and £o and not 
to the interfacial total charge change of sign, like stated in the former author's work [7J. 

It is interesting to stress that the hydraulic model has the added value of being simple and intuitive. This has 
allowed us to find, for the first time, the stabilizing behaviour of polarization charges in the ohmic/ohmic interface, 
even though this type of interface has been extensively studied (this result is also valid for the infinite interface). In 
the non-omhic/ohmic interface, the model also allowed us to find out the true reason for the appearance of stability 
bands in the non-dimensional formulation in a previous work 0. Thus, the hydraulic model allows us to find new 
unknown features and to correct mistaken interpretations in previously studied systems. The hydraulic model also sets 
a reference frame that can be used to find out for what values of the system parameters (conductivities, ion mobilities, 
dielectric constants, relative thicknesses, etc.) the electric pressure is acting as a destabilising mechanism in a two 
layer fluid interface. In addition, the model and its relation to the standard problem of a long wave perturbation in 
an infinite interface, has been set in a very formal and general approach. This allows for the model to be used as a 
first step to study a variety of interfacial (linear and non-linear) stability problems. 
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